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Abstract 

A flat elastic sheet may contain pointlike conical singularities that carry a metrical "charge" 
of Gaussian curvature. Adding such elementary defects to a sheet allows one to make many 
shapes, in a manner broadly analogous to the familiar multipole construction in electrostatics. 
However, here the underlying field theory is non-linear, and superposition of intrinsic defects is 
non-trivial as it must respect the immersion of the resulting surface in three dimensions. We 
consider a "charge-neutral" dipole composed of two conical singularities of opposite sign. Unlike 
the relatively simple electrostatic case, here there are two distinct stable minima and an infinity 
of unstable equilibria. We determine the shapes of the minima and evaluate their energies in the 
thin-sheet regime where bending dominates over stretching. Our predictions are in surprisingly 
good agreement with experiments on paper sheets. 

PACS numbers: 68.55.-a,46.32.+x,02.40.Hw 



The first step towards an understanding of a linear field theory is often to examine its 
fundamental singular solutions. These serve as elementary components of complicated fields 
constructed via linear superposition; familiar examples abound in electrostatics pQ and 
viscous [2] and inviscid [3J fluid flows. 

There is currently no analogous procedure for generating the shape of an elastic sheet from 
singularities in its Gaussian curvature. The corresponding field equations are nonlinear, and 
the canonical singularities may be more exotic than points or lines [IHH]. Yet in a thin sheet, 
the simplest defects are pointlike and conical, and such a defect bears a strong resemblance 
to a monopole. 

The intrinsic geometry of a conical monopole is characterized completely by an angle If 
we remove a wedge of material < 0), the only bending equilibrium of an isolated sheet is 
a circular ice cream cone. If this deficit is replaced by a surplus > 0), an infinite number 
of dihedrally symmetric equilibrium shapes exists [9], but the two-fold symmetric state is 
the only stable one [10J. 

The question of how to construct multipolar analogues of a conical defect is an obvious 
next step, one that has already been hinted at in the treatment of disclination dipoles as 
"dislocations" in crystalline membranes [11] . Composite defects can always be constructed 
by placing additional pointlike defects on the sheet. Indeed, when a cardboard box is folded 
from a flat template, positive Gaussian curvature will be generated at its eight vertices. 
More generally, the defect will consist of a combination of deficits and surpluses. The 
simplest composite defect of this form is a dipole consisting of two points separated by a 
finite distance, with a deficit at one point and an equal surplus at the other. Such defects 
arise in metric-preserving deformations of flat sheets by various types of forcing [T2HT6] . The 
distance, and charges, must be finite as we have already hit a snag in our construction: the 
classical dipole limit is inaccessible as we are unable to remove an arbitrarily large wedge 
from a plane. A greater complication arises from the absence of linear superposition. It is not 
generally possible to sew two cones together without introducing stretching somewhere on 
the surface. One resolution, which we will use in this work, is to approximate the geometry 
by two conical parts joined by a ridge, with planar faces on either side interpolating between 
the cones. 

We will begin by describing this geometry. The energy associated with the sheet is bending 
energy, quadratic in the mean curvature. We will ignore the energy associated with the 
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(a) (b) 

FIG. 1: Dipole geometries modeled with an A4 sheet of paper, (a) Torque dipole and (b) force 
dipole of strength ^ = ? and dipole length D = 6 cm. 

formation of the ridge between the two defects. We will see that — unlikely as this may 
appear — this is legitimate for understanding the broad features of a dipole. Even this 
simple model displays behavior without any electrostatic analogues. There are two distinct 
stable dipole states with very different energies, separated by a high energy barrier; the 
ground state is quite rigid. 

A simple way to familiarize oneself with the peculiarities of a dipole in an unstretchable 
sheet is to tinker with a sheet of paper [TT| . Take two points near the center of a rectangular 
sheet and connect them with a line parallel to the longer side. To construct a dipole of 
strength first cut out a wedge with an angle ^ at one of the two points, join the two 
exposed sides in the sheet, and round off the edges. Then, make a cut from the second point 
to the boundary, insert the wedge, and round off again. 

Two different stable shapes are possible (Fig. [I]): either the two sides extend into two opposite 
regions of space, or the whole sheet lies beneath the line joining the two points. The two 
configurations are stable with respect to small deformations, but can be transformed into 



one another if significant forces are applied. Outside the singularities, the dipole in Fig. 1(b) 



resembles the zero-deficit "developable cone" obtained, among other ways, by depressing a 



sheet into a cup with a point force at its center fT2[ ITEHiZU] . The second dipole, in Fig. l(a 
seems to represent the analogous application of a point torque. Accordingly, we will refer to 
the two configurations as force and torque dipoles. 
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In both cases, the surface can be approximated by planar parts which connect a deficit to a 
surplus cone. The two planar regions of the torque dipole appear to form triangles extending 
to the boundary of the sheet. The two triangles meet along a straight ridge, with dihedral 
angle 2a, connecting the conical singularities. The force dipole can be approximated in 
similar fashion, but one must exercise a little more care: in the paper model, the ridge is not 
straight, and it is not completely clear how the flat regions can be defined. We will return 
to this point below. 

For the theoretical model, we will assume that the dipole consists of a pair of unstretchable 
cones joined by a straight ridge and planar regions. The dipole vector D points from the 
deficit singularity to the surplus singularity along the ridge. This vector, together with 
the apparent mirror plane of the equilibrium dipole geometries, defines the perpendicular 
vertical z. 

To construct the global geometries, we first consider the two conical defects separately. We 
will use the subscript ± in the following only when the distinction between the two cones is 
relevant. Each of the cones can be described in terms of an open curve r± : s — > u±(s) on 
the unit sphere, where the arc length s runs from to L±. The value of L± depends on the 
dipole strength ^ and, implicitly via the coupling between the two cones, on the dihedral 
angle 2a. Each conical surface is given by the vector function X±(r,s) = ru±(s) G M 3 , 
where r is the radial distance to the apex. The two tangent vectors to each cone, u and 
t = d s u = u' , form a right-handed orthonormal surface basis together with the normal 
n — u x t. Whereas the surface is flat in the direction of u, its curvature along s is given 
by k = —n ■ t' . The latter corresponds to the geodesic curvature of T on the unit sphere 
and should not be mistaken for the Frenet curvature. 

We will have to identify the configurations that minimize the bending energy of an open 
unstretchable cone of radius R with surplus or deficit angle. Since tangent vectors and 
curvatures must be continuous across the boundaries between flat and conical parts, k must 
vanish along two fixed radial lines. If we cut off the apex at r$ and integrate over the 
radial direction, the bending energy is given by B = (a/2) J T dsn 2 , where a = \n(R/ro). 
The unstretchability of the sheet is implemented by adding a constraint term to the energy 
functional which fixes the metric of the surface [2T] . 

For an isolated conical monopole, there is a remarkably simple way to determine the shape 
of such a surface [HI EE]- Recall that one can find a conserved Noether current for every 
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continuous symmetry of a system. The fixed apex breaks translational invariance of the 
monopole, but the bending energy remains rotationally invariant. The corresponding con- 
served vector M is given by M = (k 2 /2 — C) n + n't + ku. Its square yields the first 
integral of the shape equation of the monopole: 

M 2 — C 2 = k' 2 + k a /A + (1 - C)k 2 , (1) 

analogous to planar Euler elastica. The constant C is associated with the stress in the 
surface due to unstretchability. 

The problem we consider here is more complicated than that of an isolated monopole: the 
surface is modeled by combining two individual cones, neither of which forms a closed surface. 
One might be skeptical as to whether it is correct to appeal to the rotational invariance 
underlying conservation of the torque vector M for each individual cone. However, we will 
see below that this yields a valid approximation for the constituent cones of the dipoles. 
The interaction between the two cones will occur entirely through boundary conditions. 
Eqn. Q can be solved for k in terms of elliptic functions. Since k vanishes on 
the boundaries, one can deduce directly that the constant of integration M 2 — C 2 
has to be positive. One obtains for the curvature (see appendix [A|: k±(s) = 
^f4n±y/—k (JC [k] / L±) sn [2n±(K, [k] /L±) s, k], where sn [s, k] is the sine of the Jacobi am- 
plitude am [s, k] with parameter k < 0, K [k] is the complete elliptic integral of the first kind 
[22], and n± denotes the number of lobes of each cone. In the paper model, n_ = 1 and 
n + = 3. The vector M± is conserved in each of the cones, but will not generally be aligned 
with the vertical axis z. Let each M be aligned along a unit vector j, with M = Mj and 



M > (see Fig. 2(a)). Due to reflection symmetry, M must lie in the mirror plane Dz. 
Let the generators along the boundary of each cone point in the directions i = u\ s=0 and 
Ri = u\ s= l ± , where R represents a reflection in the mirror plane. For each of the two cones, 
the vector i will be parameterized by two angles: the angle (3 G {0, 7r} it makes with the 
dipole vector and the angle a e {0, 7r} the Di plane makes with the mirror plane Dz (see 



Fig. 2(b)). When the dipole strength is zero, a = tt/2 and the sheet is flat. The angle a 
must be the same for both cones. Given \l/ and the ratio D/R, the geometries we consider 
are completely characterized by a, which will be chosen by minimizing the total bending 
energy of the sheet. 

The length L of the curve F depends on the angles /3 and \l/: L± = 2n — 2/3±±\]/. It is evident 
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(a) (b) 

FIG. 2: Definition of variables for the surplus cone. The same conventions are used for the deficit 
defect. 

that \T/ cannot exceed 2ir, a bound on the deficit and, thus, the strength of the "charge- 
neutral" dipole. We cannot take the dipole limit \1/ — > oo, D — > with ^>D = const. > 0. 
There is no exact analogue of an ideal electrostatic point dipole. 

As explained in appendixlBl we use the projections of M onto the local trihedron (u, t, n) to 



determine f5± as a function of a [26]. The results for the torque dipole are shown in Fig. 3(a) 
for several values of \T/. For a fixed \T/ one only finds solutions for /?_ up to a maximum value 
of a. Similarly, there is a minimum value of a below which no solution for /3+ can be found. 
For \T/ — > 0, both of these limits go to |, which corresponds to a planar sheet. Solutions for 
a > | can be found by inverting the sign of the curvature k±(s). This simply corresponds 
to a rotation of an equilibrium solution around the dipole vector D by an angle of it. Since 
such an operation does not change the angle (3, one finds (3±(tt — a) = (3±(a). This relation 
is also true for the force dipole. If we again look at a < f, the curves of /3_ are the same as 
for the torque dipole. With the original choice of signs for k + (s), one only finds solutions 
for the surplus cone of the force dipole when a > |. By rotating about D, one can obtain a 



plot corresponding to Fig. 3(a) (see Fig. 3(b) ): one again finds a minimum value of a which 



is now determined by the requirement that we only want to consider shapes which do not 
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FIG. 3: Upper panel: The angle /3 as a function of a for (a) the torque dipole and (b) the 
force dipole, with strength ^ = ^ (thick blue curves), | (green curves), it (red curves), and ^ 
(thin black curves). The dashed curves show /3_(a), the solid curves /3+(a). Shapes where the 
surplus part of the force dipole intersects with the planar regions correspond to the dotted part 
of the curves for /3+(a) in (b). For the torque dipole, intersections occur for higher values of ^ 
0& > 5-14). 

Lower panel: The shaded zone represents those values of (vP, a) for which values of /3+ and /3_, and 
thus solutions, may be found for (c) the torque dipole and (d) the force dipole. Different colors 
represent qualitatively different cases. For example, in the blue region /3_ > | and /3+ < |. The 
double curve indicates shapes where the surplus part touches (c) the dipole vector, (d) the planar 
regions of the dipole. 

7 




7T 

2 



a 



■- 

4 

mm 



- 

4 



(a) 



(b) 



FIG. 4: (a) Minimal bending energy as a function of \l/ for the torque (thick curve) and 

the force dipole (thin curve). Insets: Total bending energy B tot of each dipole as a function of a 
for $ = ?. (b) Dependence of the angle a m ; n on ^ determined semi-analytically (solid curves) 
for both dipoles, and experimentally for the torque dipole using an A4 sheet of paper with dipole 
vectors D = 2cm (red triangles), 4cm (green squares), and 6cm (blue circles). Points and error 
bars are mean and standard deviation of four measurements. 

make self-contact. In fact, several types of self-contact are possible: (i) the surplus cone 
can touch the planar regions of the dipole, (ii) adjacent folds of the surplus cone can touch 
each other, (Hi) when D/R is too small, surplus and deficit cones can touch each other in 
the force dipole configuration. Whereas we can exclude case (Hi) by assuming that D/R is 
sufficiently big, we will have to consider (i) and (ii) here. Numerically, one finds that the 



surplus part touches the planar regions before it touches itself (double curve in Figs. |3(c) 



and 3(d)). Results found for the monopole suggest that shapes with self-contact will break 
symmetry |23j . 

We determine the shape of the dipole by combining these isolated solutions. The angle a has 
to be the same for both parts. In general, the region in the (\l/, a) space in which a solution 
can be found is restricted, as can also be seen in the lower panel of Fig. [3j To determine 
which a is predicted by our theoretical model, we minimize the total bending energy of the 
dipole. We neglect the influence of the width of the sheet and the length of the dipole vector 
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D, and consider B tot = B + + B_, where B± = J L± ds4/2 = &n%K, [k] (S[k] - K [h])/L±. 
The function S[k] is the complete elliptic integral of the second kind [22J. The insets of 



Fig. 4(a) show B as a function of a for ^ = f . For both dipoles we find one minimum 



-B^jj at an angle a m i n (circular dots in the insets). Interestingly, the sum of the conserved 
torque vectors M tot = M + + M _ is vertical at this minimum. The qualitative behaviour 
of B tot does not depend on for the torque dipole, the minimum always lies somewhere 
in the middle of the curve, whereas for the force dipole, it coincides with the upper limit 
where /3_ = 0. 

These results help us to understand the behavior of the paper dipole (see Fig. [I]). We observe 
that it is simpler to convert a force dipole into a torque dipole than vice versa. In terms 
of our theory, this can be explained by the difference in energy of the two states: for fixed 



the torque dipole always has a lower energy than the force dipole (see Fig. 4(a)). Since 
both conformations are stable with respect to small deformations, there must be an energy 
barrier between them. 

Moreover, we can now compare the theoretical value of ct m i n to experiments. For the force 
dipole, our theory predicts that a will be as large as possible. However, the real physical 
sheet can also stretch outside of the singularities, and thereby increase the dihedral angle 
even further. Consequently, the real ridge is not entirely straight, as mentioned earlier. The 
angle a min seems to be approximately equal to | close to the ridge, but its value depends 
on where along the ridge it is measured. In the paper model of the torque dipole, the angle 
a m i n is approximately constant along the ridge. We constructed four paper models of the 
dipole for each of five strengths \l/ and three dipole lengths D. The angle a m i n for each 
model was obtained using one of two different methods. In one, we placed a compass on the 
ridge and subsequently measured the angle with a protractor. In the other, we determined 
tan(a min ) as given by half of the wingspan of the model divided by the vertical height of the 
ridge. The results are plotted in Fig. 4(b)[ One finds that a m - m increases with D for fixed 
a dependence that our model does not consider. This effect is, however, small compared 
to the functional dependence of a m \ n on the dipole strength Our theory explains this 
functional dependence rather well. 

We have found, via paper constructions, two stable local minima of bending energy for a 
dipole in an elastic sheet. We have assumed that each of these states can be decomposed 
into two conical singularities joined by flat parts. This apparently crude assumption explains 
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the resulting non-trivial shapes surprisingly well. The force dipole we have constructed is 
reminiscent of a well-known fundamental component of crumpling [T5], while its torque 
counterpart resembles a numerically-generated buckled dislocation in a crystalline sheet 
[TTJ . The rigidity of the latter recommends it as a potentially useful element in lightweight 
load-bearing structures. For either construction to be incorporated within a larger system of 
defects, a closer look at their far-field asymptotic geometry will be necessary. It is possible 
that other structures may be induced by the interaction of the constituent monopoles [24J . 
The role of this interaction in the mechanical response of the surface is another relatively 
unexplored topic [TJ1 [TB], although in this case the distinction between our constructions 
and the more mobile analogues induced by forcing may be important. 
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Appendix A: Determining the curvature 

If we assume C < 1 and write A -2 := 1 — C, the first integral of the shape equation 



M 2 -C 2 = k' 2 + - + {1-C)k 2 



(Al) 



can be integrated: 




where M 2 - C 2 > . 



(A2) 



Inversion yields 




where rj = \f \ + (M 2 — C 2 )A 4 > 1. The function sn [s, k] is the sine of the Jacobi amplitude 
am [s, k] with parameter k [22J. 
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The constant A can now be determined by requiring that k(0) = k(L) = 0. There is, 
however, a qualitative difference between the two defects: the curvature of the conical part 
with the deficit angle has one extremum, while the curvature of the part with the surplus 
angle has three. Consequently, one obtains for n + = 3 and n_ = 1: 



A 



± 



2n±>/2K, 



1-7? 
l+ri 



(A4) 



The symbol K. [k] denotes the complete elliptic integral of the first kind [22] . Without loss 
of generality we assume that the curvature of the deficit part has one maximum and no 
minimum whereas the curvature of the surplus part has two minima and one maximum. If 
we define k :- 



\^ e [-1,0], we obtain: 



K±(S) 



An 



-kK[k] 



sn 



2n ± /C [k] 



s, k 



(A5) 



The solutions with the opposite sign, i.e., with the opposite extrema, can be obtained by 
rotating the resulting shapes by an angle of n around the dipole vector. Since X" 2 := 1 — C, 
one has: 

4n|(l + fc)/C[fc] 2 
C± = 1 — 2 . (A6) 

A similar calculation can be done for C > 1. One finds the same expressions for n±{s) and 
C, but now k = 1±= < -1. 



Using the quadrature (Al) at s = 0, M — a/k'(0) 2 + C 2 , one can determine the absolute 
value of the torque vector M in terms of k and L±: 



M 4 



i 



64n 4 ± k!C[kf 
Li 



+ 



1 



4n|(l + k)IC[kf 



(A7) 



The only parameters left undetermined are k and L. They are given implicitly by the 
boundary conditions at i and Ri (see next section). 



Appendix B: Identifying the two conical shapes 

To find the value of for fixed angles \1/ and a, we use the projections of M onto the local 
trihedron (it, t, n) (see Fig. 2 of the main text). The direction of u in Euclidean space is fixed 
by its azimuthal and polar angles on the unit sphere, if and •& respectively. The azimuthal 
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/3<f 


/3>f 


a< § 




+ 


a> § 


+ 





C7_ 


/3<f 


/3>f 


«< f 


+ 




«> f 




+ 



TABLE I: The sign of C depends on the value of a and (3. 



<t>+ 


/3<f 


/3>f 




0+ <o 


0+>O 


0+ >0 


<t>+ <o 


«<! 


n.p. 


2ir — 7 


2-7T — 7 


self-contact 


a > § 


self-contact 


7 


7 


n.p. 





/3<f 


/3>f 




0_ < 


0_ > 


0_ > 


0_ < 


a< f 


7 — 2ir 


7 


7 


n.p. 


a> f 


n. p. 


n. f. 


n. f. 


n.f. 



TABLE II: The functional dependence of <j>± = ip±(L±) on the angle 7 depends on the sign of <j>± 
and the value of a and f3 (see text). Shapes with both C± < and <p± < are not possible (= n. 
p. in the table). For the deficit cone one does not find solutions if a > J (= n. f. in the table). 



angle is measured with respect to i and the polar angle with respect to j. Projecting 

,2 



M 



K 

~2 



C I n + n't + ku 



onto u 



M(j ■ u) = Mcostf = k 



(Bl) 



(B2) 



provides the polar angle 1? as a function of the curvature k. Eqn. (B2) places a strong 
constraint on the equilibrium shape. In particular it implies that the radial lines where k 
vanishes lie on the plane perpendicular to M, i.e., i _L j and Ri _L j. It also implies that 
the curvature k is strictly positive above this plane and strictly negative below it. 



Projecting M onto t reproduces the derivative of Eqn. (B2). The remaining projection 

1 



M(j ■ n) = Mtp' sin 2 1? = -kt 



C 



(B3) 



determines the azimuthal angle tp in terms of the integrated curvature 



1 



Msin 2 ^ V 2 
It is known that [22 



K 2 -C 



\k 2 



C 



-M 



M 2 ' 



M 2 -2C 

W 2 



1 - M~ 2 H? 



(B4) 



ds 



1 — 6 sn 2 [cs, k] c 
12 



IT [b, am [cs, k] , k) 



(B5) 



where II is the elliptic integral of the third kind. Thus, with Eqs. (B4) and (B5) 
V?±0) = / dt(p'(t) 



Ml - 2C 4 



2M± 2n±K [k] 



n 



16n|( 



-k)lC[kf 



, am 



2n±K, [k] s 



k 



■>± 



k 



M 4 



(B6) 



s. (B7) 



Introduce a Cartesian chart on the surplus cone, adapted to the geometry of the dipole. 
The unit vector x is parallel to the dipole vector D and points away from the ridge (see 
Fig. 2 of the main text). This, together with a vector y and the vertical vector z, defines a 
right-handed basis. Likewise, one can define a different Cartesian chart on the deficit cone 
with a unit vector x anti-parallel to D. 
At s = the trihedron (it, t, n) is then given by: 



i — — cos (3 x + sin (3 i± , t s= o = sin /3 x + cos /3i± , n 



s=0 



— cos ay + sinaz , (B8) 



where i± = — sin ay — cos a z. The unit vector j in the direction of M can also be expressed 
in terms of the Cartesian basis (x, y, z): 

i x Ri cos a sin x — cos z (B9) 



I* x R*| v 7 ! - sin 2 a sin 2 /3 



The projection of M = Mj onto t has to be equal to k' ± (see Eqn. (Bl )). With the definition 



of signs in Eq. (A5), one obtains the sign of j as follows: 

.«<5 i x Ri 



J± 



, and 2 



J± 



\i x Ri 

since we have defined M > 0. The projection at s = thus gives us for all a G {0, tt}: 



(BIO) 



M-t |cosa| 



'1 — sin 2 a sin 2 (3± 
Projecting M onto n at s = yields: 

M± sin a cos j3± 



c a< f 



3 



and (7° 2 



(Bll) 



(B12) 



a/1 — sin 2 a sin 2 (3± 

The sign of C depends on the value of a and (3 (see Tab. ffl): signC-t = 
Tsign [(f-«) (f-/5)]. 

Finally, we have to connect (p(s) to the extrinsic geometry. The range of the azimuthal angle 
= <f{L) on the iRi plane depends on the angle 7 between i and Ri in a nontrivial way 
(see Tab. [TT|). 7 is given by 



cos 7 = i ■ Ri = cos 2 (3 + cos (2a) sin 2 (3 , where 7 G {0, 7r} , 



(B13) 
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whereas one finds the following expression for with the help of Eqn. (B7): 

-k)K [Jfef 



2C 4 



n 



M 2 ± L\ 



k 



M 4 



(B14) 



2/C [k] M± 

where we have used that am [2n±K. [k] ,k] — n±7t for k < 1. The value of can in principle 
be positive or negative. For > 0, one either has = 7 if j and i x Ri are parallel, or 



: 2tt — 7 if the two vectors are anti-parallel (see Eqn. (BIO)). When < 0, one finds 
7 — 27r if j and i x Ri are parallel, and = —7 if they are anti-parallel. Using Eqn. (Bll ) 



and the condition on (see Eq. (B14) and Tab. [TTJ) one can now calculate (3± as a function 
of a. 

It is instructive to discuss the sign of in connection with the sign of C (see Tab. [T]). If C 



is negative, Eqn. (B4) directly tells us that <p'(s) > for all s E {0,L ± } and consequently 
> in this case. This implies that both C < and < are not possible (= n.p. 
in Tab. [Tl]). For C positive, the situation is more complicated. <p' at s = is now always 
smaller than zero, but changes sign along r if C is not too big. To determine C cr , where 
changes sign we can, for instance, consider the projections of M at s = L/2 where the 



curvature k is maximal: <f>± changes sign when d(L/2) = 0, i.e., when M° 



k(L/2) and 



2C C 



k 2 (L/2) = (M 



cr\2 



Appendix C: Harmonic approximation 

Let us suppose that ^//2n 1, and expand the curvature k(s), as well as the constants C 
and M, in powers of e = \/—k: 



k = Ki + k 3 + 



M = M + M 2 + ■ ■ ■ ; C = C + C 2 + 



(CI) 



For \J— k — » 0, one has /C [0] = | and sn [m s, 0] = sin (m s). The equilibrium states are de- 
scribed by small oscillations about an equator, where k = 0. In the harmonic approximation 



of the quadrature (Al) about k = 



M 2 - C 2 + 2M M 2 - 2C C 2 = (k',) 2 + (1 - CoK • 



(C2) 



At lowest order, the arc length coincides with the azimuthal angle, if. The curvature along 
a loop, vanishing at ip = and at L , is given by 



\Ki\ & Ai sin I 



(C3) 
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where the amplitude A\ is a constant. From Eq. (A5) we infer: A\ 



2nn . / E 
L V «• 



The quadrature implies the two zeroth order identities 



M 2 = C 2 , 



(C4) 



where M is positive, as one can show by expanding Eq. (A7). The sign of Cq depends on 
whether nir is bigger or smaller than L (see below). 



Eq. ( C2 ) also implies 



2(M M 2 -C C 2 ) = Al(l-C ) 



and thus 



-2C n 



-M 2 sign (C ) + C 2 ] . 



(C5) 



(C6) 



1-<V 

This is a constraint on a combination of the second order corrections M 2 and C 2 . As a check 



of consistency, one can expand the nonlinear expressions for C and M, Eqs. (A6) and (A7), 



in v— k and show that they indeed obey Eqs. (C4) and (C6). 

The amplitude A\ depends on the ratio of arc length to polar angle over the trajectory on 
the unit sphere. One has the exact formula 



L ± = 2tt - 2(3 ± ^ 



where \<f>\ = 2n — 7 if (3 < | or 



dip 



— + sin 2 1? 
dcp J 



1/2 



(C7) 



7 if (3 > f. Using Eq. (B2) and the harmonic 



approximation for k given by Eq. (C3), one obtains 



sin 2 $ « 


1 - 




1 


dp} J 


Ml 



Ml 



sin 2 — C , and 



:i-gp) ^ 

M 2 



cos 



l-C p 



(C8) 
(C9) 



Thus, in the quadratic approximation, 

2tt - 2/3 ± ^ « |<; 

Note that \</>\ = L + <D(\f-k). Together Eqs. (|C6j) and (IClOj) imply the identity 



1 - 



4C 



(CIO) 



- M 2 sign C + C 2 = 2(1 - C )(r - 1) 
where T = (2tt - 2/3 ± $)/\(f>\. 



(Cll) 
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This is all we need to satisfy the boundary conditions at lowest non-trivial order [25J. 



The remaining boundary condition (|B12|) only involves C and M in the combination 

C 2 -M 2 sign (C ) 



-M 2 sign (C ) + C 2 at second order, or 



C C M C 2 - C M 2 . 
M^Mo + =Slgn(Co) 



M 2 



1 + 



C 



Using Eq. (Cll) this reads 

C_ 
M 



sign (C 



i + ^(r-i; 



(C12) 



(C13) 



We now approximate the right-hand side. To this end, we need an expression for the angle 
7. When is small, a ~ 7r/2, and it is possible to approximate 

2 



cos 7 « cos (2/3) + 2 ^- - a J sin 2 /3 



or, equivalent ly, 



7 = 2/3 -A7. 



if /3 < 



7 = 2vr-2/3 + A7, if /3 > 



7T 

2 

7T 



and 



(C14) 

(C15) 
(C16) 



where we have defined A7 = (| — a) tan/3. 

For both /3 < | and /3 > |, we thus have |0| = 27r — 2/3 + A7 and, 

1 



T - 1 



2tt-2/3 



[±* - A7] 



Furthermore, using Eq. (C4) we obtain 

2 



C, 







l-C 



1 







2tt - 2/3 



nn 



2A 7 



2tt - 2/3 
n 2 7r 2 



(C17) 



(Cl{ 



the correction in 7 is irrelevant at this order. One thus has 



2(1 -C ) 

Cn 



(r-i) 



2tt - 2/3 



2tt - 2j9 



(=F* + A 7 ) 



(C19) 



When a ~ |, the trigonometric function on the right-hand side of Eq. (B12) reduces to 



sin a cos /3 / n 

^^^^^^^^= « sign p 

1 -2 ■ 2 a V 9 

1 — sin a sin /3 



(C20) 



To formulate the boundary condition (B12) in the correct manner, we have to take a closer 
look and determine the signs of the right-hand and the left-hand side of the equation for the 
different cases. 



16 



We know that Lq has to be smaller than 2ir. Since n + = 3 for the surplus part, we directly 
obtain from Eq. (C4) that C +i o < — § is always negative and M = — C +t o. The (nonlinear) 
C + is thus also negative for small one only finds two types of solution (see Tab.[lJ: (a < |, 

/3 + < f , L = 2tt - 7) and (a > § , /3+ > f , L = 7). 



Consequently, the boundary condition (|B12|) reads 
2tt - 2/3 



2tt - 2/3 



1 - 



3tt 



a I tan (3 



1 /7T 

2 V2 



a) sec 2 /3. (C21) 



for all /3+ G {0,vr}. 

For the deficit part n_ = 1. Since for /3_ < | we always have L = 2rr — 7 > tt, we find 
that C_ )0 is positive. For /3_ > | we always have L = 7 < 7r and consequently C_ >0 < 0. 
Both types of solution with a < | are thus allowed (see Tab. [Ijfor the sign of C_). However, 



a > I is not allowed. 



The boundary condition (B12) for /3_ G {0,7r} becomes: 

,n -1 



2tt-2/3 



2tt -2/3 



7T 



— — a J tan /3 



) =~\{l -«) 2 sec 2 /3. (C22) 



Eqs. (C21) and (C22) can be reorganized in the following form: 

1 „ /7T 



where 



/±C9) 



±- (2tt - 2/3) sec 2 /3 



±- sec 2 /3 



(2tt - 2/3) 



a: 



2tt - 2/3 
n±7r 

2tt - 2/3 

7l±7T 



(C23) 



± 2tan/3 
+ 2 sin 2/3 



(C24) 



For a nonvanishing Eq. (C23) implies that a = f is only possible for values of /3-t at 



which f± diverges. The function / + of the surplus part diverges at /3 = | whereas the one 
for the deficit part is a monotonically decreasing function of (3 with no divergences. Thus, 
(3\ a= K = f for any whereas P\ a= K = only if ^ — >■ 0. The same behavior is found in the 
nonlinear regime. Moreover, the harmonic approximation and the numerical results of the 
full nonlinear equations coincide in the limit of vanishing 
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The angle a has to be the same for the two parts. For fixed \l//27r <C 1, this requirement, 
together with Eq. (C23), implies that / + (/3 + ) = in the harmonic approximation. 
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